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Abstract. We compare a multi-level statistical model with a two-level model for the decay out of superdeformed rotational 
bands in atomic nuclei. We conclude that while the models depend on different dimensionless combinations of the input 
parameters and differ in certain limits, they essentially agree in the cases where experimental data is currently available. The 
implications of this conclusion are discussed. 



INTRODUCTION 

Superdeformed (SD) rotational bands occur in the second minimum of the potential energy surface in deformation 
space. They have been observed in several mass regions around A=2Q, 40, 80, 130, 150, 165, 190 and 240 |1]. 
Characteristic of the decay of most superdeformed bands is the sudden disappearance of the total intra-band decay 
intensity when a certain spin is reached. The intensity reappears as electromagnetic transitions in the normal deformed 
(ND) minimum. In certain special cases the decay path subsequent to the SD band is known completely as it is 
for 133 Nd (3, or almost completely as for 59 Cu |3]. More typical however is the situation in the A = 80, 150 and 
190 regions, where the decay from superdeformed to normal states is spread over many different available paths, 
making observation of discrete y rays linking SD and ND states very difficult. Because of this fragmentation of the 
SD intensity, experimentalists have only been able to identify a small number of the strongest paths in a few cases 

I5I Icj I7I la. I9I fiol [Till and these only when very large data sets have been obtained. In these more typical regions 
it is reasonably clear that there remains barrier between the minima at the decay-out spin and consequently that the 
decay-out occurs by tunneling HQ Ed. 

The fact that the decay path is unknown has been thought to indicate that a statistical treatment of the ND states and 
their coupling to the SD band is appropriate. This approach is expected be valid when the ND level density is 

so high as to make a useful description of the coupling to individual ND states unfeasible. However, given that there 
is uncertainty concerning the ND level density around the SD band there remains the possibility that the decay-out 
occurs entirely through a single ND level. The ND level density depends strongly on the excitation energy of the SD 
state; as SD bands have been observed at relatively low energies in the Pb isotopes, the point of view encompassed by a 
two-level model where the decaying SD state couples to a single ND state (itself decaying) could be more appropriate. 
Indeed, in Ref. 1 19] the decay-out in the A=190 region is attributed to the crossing of the SD band with the nearest 
neighboring excited ND band. 

In this contribution we propose to compare the two limiting descriptions mentioned in the previous paragraph 
as embodied by the multi-level statistical model of Refs. 1201 12III which exploit analogies of the SD decay with 
the theory of compound nucleus reactions [ 22] and the two-level model of Refs. l23l l24l l25Tl . Neither approach 
provides a microscopic description of the decay-out such those provided by the cranked Nilsson-Strutinsky model 
of Refs. 1 12ul3lll4l.ll5ll . the cluster model of Ref. 1 19] and the generator coordinate model of Ref. 12611 . Rather, both 
models considered here are phenomenological models which attempt to provide simple formulas for the total intra- 
band decay intensity in terms of the most relevant parameters. Instead of analyzing the experimental data, in what 
follows we shall only compare the statistical and two-level models in their formal aspects. 
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FIGURE 1. Schematic diagram of the multi-level statistical model. The decaying SD state under consideration has energy £$ 
and spin J and mixes via the coupling Vjs with ND states of the same spin whose energies, Ej, are the eigenvalues of the GOE 
Hamiltonian if aw' - The SD (ND) states have a common electromagnetic width T$ (TV). 



MULTI-LEVEL STATISTICAL MODELS 



Let us denote the decaying SD state by \S), its energy by s$ and its electromagnetic width by F s . It is coupled to 
the ND states |7V), N — 1, ...,K by Vns- The ND states are modeled by the GOE Hamiltonian H NN i and are assumed 
additionally to have a common electromagnetic width Fn- The intra-band decay (see Fig. [3 is described by the Green's 
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The total intra-band decay intensity is then given by 
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and the total decay-out intensity by 
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where P{H NN i) is the probability distribution of the GOE and d[H NN i] is the product of the differentials of all 
independent matrix elements of H NN i 12^1 . Gu and Weidenmuller calculated F$ 1 20] by adapting the result for the 
ensemble average of the compound nucleus cross-section derived by Verbaarschot et al. 12911 to obtain, in the limit 
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The density of ND levels is denoted by D. From Eq. Gu and Weidenmuller deduced the fit formula 



F n - 
r S — 



1- 0.9139 (F N /D) 



0.2172 



exp ■ 



0.4343 In (T /F s ) - 0.45 (F N /£>) 



-0.1303 



-0.1477 



(7) 



(8) 



(9) 



(10) 



Energy averaging provides an alternative approach to calculating the average values of observables which depend 
on a random Hamiltonian. The energy average of G(E) is given by 



G(E) = / G{E')p{E,E')dE' 



(11) 

where the smoothing function p(E,E') is normally a Lorentzian or a box function |22]. The energy average is carried 
out for a single realization of the GOE Hamiltonian H NN t and is expected to be equal to the ensemble average to the 
extent that the GOE is ergodic. In practice calculations proceed by choosing a representation - the optical background 
representation of Kawai, Kerman and McVoy 113 Oil - which is defined such that the couplin gs V ns have statistical 
properties which are convenient for analytical calculation. (It was checked numerically in Ref. 13111 that the properties 
of the Vns which are normally assumed do indeed obtain from the underlying random Hamiltonian.) In order to obtain 
analytical results it is further necessary to assume that Fn/D ^> 1 which is the principal limitation of the energy 
averaging technique. The optical background representation was used in Ref. ll2lll to calculate F$. The result has the 
same form as Eq. (J6J with F£ v still given by Eq. Q but now 
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An advantage of the energy averaging technique is that it also yields an analytical expression for the variance of the 
decay intensity Pill : 
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From Eqs. 0, 0, (jSJl and dl2l it is seen that F$ depends only on the two dimensionless variables F/F$ and F^/D. 
However it is possible to construct three independent dimensionless variables from the input variables F, F$, F^ and 
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FIGURE 2. Schematic diagram of the two-level mixing model. The decaying SD state under consideration mixes via the coupling 
V with a single ND state of the same spin and energy £n = £$ + A. 



D. In Ref. l2(ill it is conjectured using supersymmetry arguments that the entire probability distribution of 

Fs depends only on T/Fs and Fn/D. In contradiction to this conjecture the expression for the variance of F$ derived 
using the energy averaging technique, Eq. dl3t . depends on the additional dimensionless variable (r^ + F) /Fn- 

At present the energy dependence of \G$s(E)\ 2 and \Gns{E)\ 2 cannot be resolved experimentally and the focus of 
theory has been on the analysis of Fs or Ffj. Normally, F$, Fn and D are estimated theoretically and a value of V is 
extracted from the experimental F$ 133 . An intrinsic problem of the statistical model in its application to the decay 
out of SD bands is that the variance of F$ is large when Fn/D <C 1 reflecting the fact that in this limit the positions of 
the individual ND states relative to the SD state are important [20]. Since Fn/D <C 1 for the cases which have been 
studied experimentally the preceding argument suggests that the error in the extracted V must also be large. However, 
from Eqs. (1131 and (II 41 it can be seem that the variance of F$ is suppressed not only for Fn/D ^> 1 but also for 
(r + Fs) /Fn 3> 1 . This is physically plausible since in the decay of an ND state to the SD state the larger the value of 
r + Fs the less important the exact position of the SD state. 



TWO-LEVEL MIXING MODEL 

In the two-level mixing model the decaying SD state of energy Es and electromagnetic width Fs is coupled to single 
ND state of energy En and electromagnetic width Fn by the coupling matrix element V (see Fig. |2j. In analogy to 
Eqs. and 0, the intra-band decay is described by the Green's function 

Gss(E) = - 2 (15) 

E-Es + iTs/2 y—r^ 

and the decay-out by the Green's function 
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The decay-out intensity, Fn, is given by Eq. © (instead of a sum over there is just one term) and an exact calculation 
yields Efll 
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From Eq. d!7l > it may be seen that Fn depends on three independent dimensionless variables which may be taken to 
be V/A, F/A and Fs/Fn- The calculations of Ref. 11911 suggest that the value of A is rather important and in fact 
determines the decay-out spin in the sense that the decay-out occurs very near to where the SD band crosses the 
nearest excited ND band. The behavior of Fs as a function of A for a model related to the two-level model is discussed 
in Ref. (3- I n tne event that A could be calculated it would be favorable to use Eq. d!7t to extract V from the 
experimental Fn- 

An alternative to calculating A is to assume that it obeys some probability distribution. In Ref. 1 25] it was assumed 
that A has the distribution 

P(A) = f ds&> s (A)P{s), (21) 







where 

P(s) = I^T™ 2 / 4 and <J»,(A) = Ie(i-M). (22) 



In Ref. 12511 it was argued that the average value (|A|) — D/4 was the appropriate value to use in the calculation of Fn- 
Alternatively, one could calculate the average of Fn with respect to the distribution P(A): 

(F N ) = j^F N (A)P(A). (23) 

Since A = — Es, Eq. d23t corresponds to an average over the eigenvalues, £,-, of the GOE Hamiltonian of the 
statistical model. Thus, while Eq. d23i is related to the ensemble average, Eq. 0, it is not identical with it since 
Eq. (|5} is an average over the H N n>, ie. an average over both eigenvalues and eigenvectors. It is clear from Eq. (123 i 
that (Fn) depends on the dimensionless variables V/D, F/D and Fs/Fn- Both Fn(A = D/4) and (Fn) are compared 
with Fn of the statistical model in the next section. 



NUMERICAL COMPARISON OF THE MODELS 

The first thing which may be noted on comparing the statistical model and the two-level model is that they depend 
on different dimensionless combinations of the input parameters V, D, Fs and Fn- The statistical model depends on 
F/F s = 2n(V /D) 2 (D/F S ) and F N /D while the two-level model depends on V/D, F/D = (F N /D)(1 + F s /F N )/2 and 
Ts /r#. In addition, it may be deduced from Eq. (I17> that Fn in the two-level model increases monotonically from zero 
to Fn/(Fs + Fn) as V is varied from zero to infinity whereas in the statistical model Fn increases monotonically from 
zero to 1 . 

A second difference occurs in the (physically not realized) limit Fn /D ~ * °°. In this limit Fn in the statistical model 
becomes independent of Fn, being given by [1 + Ts/r] -1 , cf. Eq. @. The two-level model in the limit that F — > °° 
yields Fn = Fn/{Fn + Ts)[l + FnFs/4V 2 )]~ 1 which clearly does not become independent of Tat as Fn /D — > 00. 

It is useful to note that when a uniform distribution for A is assumed then 
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FIGURE 3. The figure shows 1-F S (thick solid), (F N ) (thin solid), (F N )v (dotted) and F N (A = D/4) (dashed) versus log 10 (V/Z>), 
calculated using Eqs. J6j, 1231 . 1241 and <17> respectively. The left column has Tn/Fs = 1000 and the right column has Iai/Ej = 1. 

The first, second and third rows have F/D = 10~ 4 , 0.1 and 10. ForF s fl we used Eq. (lOj when F/D = 10~ 4 , 0.1 and Eq. |l2j when 
F/D = 10. The negative values of Fn in the first row are due to the fact that Eq. 1101 is an approximation to the exact result, Eq. JS}. 
Eq. 1101 also appears to overestimate F$ in the middle row. Note that (Fn) and Fn(A = D/4) coincide in the bottom row. 



which reduces to F/F$ in the limit that V <C Tn,Fs. The limit of small V was considered for the statistical model in 
Ref. 1 20] where a perturbation expansion of Eq. <2i was carried out and the ensemble average of Fn calculated: the 
result is also F/F$. 

In Fig.[3]where we plot Fn versus V /D, the considerations of the preceding paragraph are born out. The graphs in the 
left hand column were calculated for Fn — lOOOTs and those in the right hand column for Fn — F$. The first, second 
and third rows were calculated for r/D=10~ 4 , 0.1 and 10 respectively. The dependence on F^/Fs is clearly seen in 
the two-level model as is its absence in the statistical model. Eq. i2Al is seen to agree with the statistical model for 
small V but to be unphysical for large V. It can also be seen that (Fn) and Fn(A = D/4) agree best with the statistical 
model when F/D is small. Indeed, when Fn 3> F$ so that the coherence effects present in the two-level model are 
suppressed andF <C D then (Fn) and Fn(A = D/4) are close to the statistical model for all V/D. It is interesting that 
(Fjv)tj agrees better with statistical model than (Fn) when f> 1 suggesting that in this region the average over the 
uniform distribution, Eq. i24\ . corresponds more closely to the ensemble average, Eq. (|5}, than Eq. (I23> does. The 
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FIGURE 4. Decay-out intensity for the three-level model. The dotted and long dashed curves are the branching rati os fo r decay 
to the nearest ND state and the next nearest ND state respectively while the thin solid curve is the their sum (see Eq. <25> and the 
discussion which follows it). The thick solid and dashed curves are the same as in Fig.|5|(they are the decay-out intensities for the 
multi-level statistical and two-level models respectively) as are values of T/D and Tn/Ts used to plot the curves. 



case T/D = 10~ 4 , Ffj/Fg = 1000 shown in Fig.[3]is representative of the mass 190 region. The considerations in the 
present contribution thus go some way to explaining why the values of V extracted [33] from experimental data in the 
mass 190 region using the statistical model and those extracted using the two-level model were very similar. 

In Ref. 1 25] the validity of the two-level model was investigated by calculating the decay-out intensity for a three- 
level model where the decaying SD state is allowed to decay to the nearest neighbor and next-nearest-neighbor ND 
states. The decay-out intensity may be written 

4 3) =£ f j> F i = ^Ej E]Gjs{E)][2 ' (25) 

where Gjs denotes the Green's function of Eq. (0 in the basis of eigenvectors, \j), of H NN i. Clearly, Fjp depends on 
A] and A2, the relative energies of the two ND states and on V\ and V2, their interactions with the SD state. In Ref. 1 25] 
a probability distribution for A2 similar to Eq. J2 1 1 was assumed resulting in the average value (A2} = 3D J A, 



In Fig. @] we compare Fn of the multi-level statistical model and Fjv(A = D/4) of the two-level model with F N , 
F\ and F2 for Ai — D/4, A2 = 3D/4 and V\ = V% = V. In similarity to two-level case, the effects of coherence are 
pronounced in F N when T^r = Tg. It can be seen that F N moves in the direction of Fn relative to Fn(A — D/4). (As 
mentioned in the caption to Fig. [5] Fn appears to be underestimated in the middle row due to use of the approximate 
Eq. dlO>.) We expect that F^ calculated analogously to F^ would rapidly approach Fn with increasing n. The fact 

that Ffp and Fn(A — D/4) are barely distinguishable for the case T/D = 10~ 4 , Tn/^s — 1000 again suggests that the 
two-level model is adequate to describe the SD decay in the mass 190 region. It would be interesting to calculate the 

average of F N ^ over Ai and A2 in analogy to Eqs. d23i and (124-i . 

CONCLUSIONS 

We have compared a multi-level statistical model and a two-level model of the decay out of a superdeformed band 
with each other. We conclude that while the models depend on different dimensionless combinations of the input 
parameters and yield different results for certain limiting cases, they essentially agree in the cases where experimental 
data are available. However, this conclusion does not mean that we have reached a fully satisfactory description of the 
decay-out of superdeformed bands. The two-level model and multi-level statistical model level quite naturally only 
agree in the parameter regime where there are just two well defined levels interacting. However, this agreement, as it is 
described in the preceding section, surely has much to do with the use of the GOE probability distribution to calculate 
the ensemble average in Eq. (0 and the use of the Wigner nearest neighbor spacing distribution in obtaining Eq. (12 1 i 
for the probability distribution of A. The use of these two intimately related distributions to account for the unknown 
parameters of the respective formulations means that level fluctuations must be equally important in the two-level and 
multilevel statistical models. Consequently, the respective values of V which can be extracted from the experimental 
Fn using the two models appear to be subject to similar errors. An analytical calculation of the variance of Fn, based 
on either model, which is valid in the physically realized region would be useful in assessing the errors. 
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